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O . Abstract 

We study the discovery potential of the Tevatron for a Z' gauge boson. We 
introduce a parametrization of the Z' signal which provides a convenient bridge 
^p Hi between collider searches and specific Z' models. The cross section for pp — ► Z'X — > 

Oh, £ + £~X depends primarily on the Z' mass and the Z' decay branching fraction into 

rH \ leptons times the average square coupling to up and down quarks. If the quark and 

lepton masses are generated as in the standard model, then the Z' bosons accessible 
at the Tevatron must couple to fermions proportionally to a linear combination of 
baryon and lepton numbers in order to avoid the limits on Z — Z 1 mixing. More 
generally, we present several families of U(l) extensions of the standard model that 
include as special cases many of the Z' models discussed in the literature. Typically, 
the CDF and DO experiments are expected to probe Z'-fermion couplings down 
to 0.1 for Z' masses in the 500-800 GeV range, which in various models would 
substantially improve the limits set by the LEP experiments. 



1 Introduction 



An important question in particle physics today is whether there are any new gauge 
bosons beyond the ones associated with the SU(3)c x SU{2) W x U(l)y gauge group. 
This question is interesting by itself, given that the selection of the gauge bosons observed 
so far remains mysterious. Furthermore, new gauge bosons are predicted within many 
theories beyond the Standard Model (SM) which have been developed to provide answers 
to its many open questions. 

The simplest way of extending the SM gauge structure is to include a second U{\) 
group. The associated gauge boson, usually labeled Z', is an electrically-neutral spin-1 
particle. If the new gauge coupling is not much smaller than unity, then the U(l) group 
must be spontaneously broken at a scale larger than the electroweak scale in order to 
account for the nonobservation of the Z' boson at LEP and run I of the Tevatron. In 
this article, we study the Z' discovery potential of the run II of the Tevatron, the highest 
energy hadron machine operating for the next few years. 

The theoretical framework for studying Z' production at hadron colliders has been 
developed more than two decades ago pQ. Nevertheless, various pieces of information 
collected recently have an impact on our attempt of addressing a number of specific 
questions: What Z' parameters are relevant for Tevatron searches? What regions of the 
parameter space are not ruled out by the LEP experiments, and would allow a Z' discovery 
at the Tevatron? In case of a discovery, how can one differentiate between the models 
that may accommodate a Z' boson? 

It is often assumed that the Z' couplings have certain values motivated by some 
narrow theoretical assumptions, allowing for the derivation of a Z' mass bound (21 El- 
The opposite approach of leaving the couplings arbitrary 4 a suffers from the existence of 
too many free parameters. However, a few theoretical constraints are sufficiently generic 
so that it is reasonable to focus on the region of the parameter space that satisfies them. 
This observation, used to define the so-called nonexotic Z' bosons underscores the 
importance of the Z' couplings to the SM fermions for collider phenomenology 6J, while 
reducing the set of Z' parameters. 

In this article, we address Z' models both from a theoretical perspective and with 
respect to their potential observation at hadron colliders. In Section El we present the 
theoretical framework needed to describe a new neutral gauge boson. We analyze the 
constraints due to gauge anomaly cancellation and the gauge invariance of the quark 
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and lepton Yukawa couplings, and discuss what new physics would soften these con- 
straints. We identify several interesting families of Z' models, and then derive the LEP 
limits. Section El is concerned with Z' production at hadron colliders, including a survey 
of theoretical tools to describe Z' events, and a convenient parameterization of limits 
from searches that simplifies comparison of experimental results with theoretical models. 
Section 4 summarizes our conclusions. 

2 Parameters Describing New Neutral Gauge Bosons 

Any new gauge boson is characterized by a mass and a number of coupling constants. All 
these parameters appear in the Lagrangian, which is constrained by gauge and Lorentz 
invariance. In this section we present a theoretical framework that is sufficiently general 
to account for the parameters that are relevant for Z' searches at the Tevatron. We 
discuss the theoretical constraints within realistic extensions of the SM, and then discuss 
the LEP limits. 

2.1 Z' mass and Z — Z' mixing 

Consider the SM gauge symmetry extended by one Abelian gauge group, SU(3)c x 
SU{2)w x U{1)y x U{\) z . The scalar sector responsible for the spontaneous breaking 
of the gauge symmetry down to SU(3)c x U{l) em includes at least one Higgs doublet and 
an SU(2)w singlet, <ft, with VEV v^. As we will explicitly show below, the constraints 
on the interactions of the U(l) z gauge boson with quarks and leptons are relaxed in the 
presence of two Higgs doublets, Hi and H 2 , with aligned VEVs vn x and vh 2 - To be gen- 
eral, we will concentrate on this case in what follows. The hypercharges of Hi, H 2 , and 
4> are given by +1, +1 and 0, respectively, so that electric charge is conserved. 

In a basis where the three electrically neutral gauge bosons, W 3tJ- , By and B%, have 
diagonal kinetic terms, their mass terms are given by: 

2 2 2 

V -f (gW^ - g Y B Y - z Hl g z Brf+ V -f (gW* - g Y B Y - z H2 g z B») 2 + ^ {z^f , 

(2.1) 

where g,gy-,g z are the SU{2) W x U(1)y x U{1) z gauge couplings, and the weak mixing 
angle is given as usual by tan^ = gy jg. The diagonalization of these mass terms yields 
the three physical states, the photon (labeled by v4 M ), the observed Z boson, and the 
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hypothetical Z' boson: 



= sin 9 W + cos 9 W , 

= W 3 » cos9 w - B» sin 9 w + eB 



= B% -e(W 3 >* cos 9 W -B£ sin 9 W ) . (2.2) 

To obtain this result we have ignored terms of order e 2 , where e is the mixing angle 
between the SM Z boson and B%, 



5M ZZ , 
C ~ Mf, - Mf ' 

The mass-squared parameters introduced here are related to the VEVs by 

,2 



(2.3) 



4 cos 2 9 



5Mh, = -j^- (z H y Hl + z H2 vl 2 ) . (2.4) 

Mz and Mz' are the physical masses of the neutral gauge bosons up to corrections of 
order e 2 . 

A Z' boson which mixes with the SM Z distorts its properties, such as couplings to 
fermions and mass relative to electroweak inputs. Precision measurements of observables, 
mostly on the Z pole at LEP I and SLC, have verified the SM Z properties at or below 
the per mil level 1 j7j , imposing a severe upper bound jH] on the mixing angle between the 
Z and Z': |e| < 1CT 3 . Therefore, it is justified to treat the mixing as a perturbation as 
was done above. 

From Eqs. ()2.3jl and ()2.4)1 it follows that the mixing angle is given by 

gz ( cos 9 W \ \z Hl + z H2 tan 2 p\ 
|C| ~ g \M 2 Z ,/M 2 Z -l) 1 + tan 2 /5 1 ] 

where tan (3 = vh 2 / v Hi ■ At least one of vh ± or vh 2 has to be of the order of the electroweak 
scale (to generate Mw, Mz and m t appropriately), so without loss of generality we can 

lr The notable exception is sin 2 9w from hadronic Z decays, which deviates at the few a level [J], and 
plays an interesting role in the fit to the SM Higgs mass 9 . These deviations have been argued as 
evidence for the existence of a Z' with non- universal interactions we shall not pursue this line of 
reasoning here. 
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set vh 2 ~ 0(246) GeV. Therefore, tan/3 > 0(1). Normalizing the largest quark U{l) z 
charges to be of order unity, the Z' production at the Tevatron is sizable only if the gauge 
coupling g z is not much smaller than unity. The mass range typically interesting at the 
Tevatron is roughly 0.2 TeV < M^/ < 0.7 TeV. Based on these considerations we find 
that the order of magnitude of the mixing angle is given by e ~ (z Hl cot 2 j3 + zh 2 ) M% / Mf, . 
The constraint |e| < 10~ 3 implies Zn x cot 2 /? + zh 2 <C 1. Although tan 2 /? could be close 
to —zh 1 /zh 2 , this would be a fine tuning, because the value of tan/? is set by the Higgs 
masses and self-interactions, and has no reason to be related to the ratio of Higgs charges. 
Therefore, in the absence of fine-tuning, a Z' accessible at the Tevatron requires |z# 2 | <C 1 
and either \zh ± \ <C 1 or tan/3 ^> 1. It is usually expected that the charges of various fields 
are either all of the same order or vanish. Although exceptions exist, such as extra 
dimensional models with brane kinetic terms jTTj on the Higgs brane, which motivate 
much a smaller effective charge for the Higgs than for the (bulk) fermions, we restrict 
attention here to the following two cases: 

z H2 = and zh x = (2.6) 

or 

z H2 = and tan (3 > 10 . (2.7) 
2.2 Couplings to fermions 

The renormalizable interactions of the Z' boson with the SM fermions are described by 
the following terms in the Lagrangian density: 

J2zf9 z Z'jYf, (2.8) 
/ 

where / = e R , tf L , u 3 R , d 3 R , q J L are the usual lepton and quark fields in the weak eigenstate 
basis; l 3 L = (v J L ,e 3 L ) and q J L = (u J L ,d J L ) are the SU{2) W doublet fermions. The index j 
labels the three fermion generations. Altogether there are 15 fermion charges, Zf. 

The observed quark and lepton masses and mixings restrict these fermion charges, so 
that certain gauge and Lorentz invariant terms can appear in the Lagrangian. In the SM, 
the terms responsible for the charged fermion masses are 

\%qid k R H + X u jk qiu k R ia^ + \ e J L e k R H + h.c. , (2.9) 
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where j, k = 1, 2, 3 label the fermion generations, and X^ k , A" fc , A^ fc are Yukawa couplings. 
In the two Higgs-doublet model described above, the Higgs doublet H is replaced by linear 
combinations of Hi and Hi- 

As discussed in section 2.1, the Z' bosons relevant for Tevatron searches have small 
mixing with the Z boson, which effectively implies that any Higgs doublet with a VEV 
of order the electroweak scale is neutral under the U(l) z symmetry. In particular, if only 
one Higgs doublet is present, then its U(l) z charge would have to vanish. Given that 
the total charge of the quark mass terms shown in Eq. (J2.9|) has to be zero, the quark 
masses and CKM elements may then be accommodated only if the quarks have generation 
independent U{l) z charges, and z u = Zd = z q , where z u and Zd are the right-handed up- 
and down-type quark charges, and z q is the left-handed quark doublet charge. One may 
relax this condition in the two Higgs-doublet model if, for example, H2 couples to the 
up-type quarks, while Hi couples to the down-type quarks, has nonzero charge, and tan f3 
is large. In that case z d may be different from z u and z q , but one still needs to impose 

z u = z q , (2.10) 

so that the large top-quark mass may be generated. 

We emphasize that we have derived this strong conclusion based on reasonable but not 
infallible arguments. One loophole is that some of the terms in Eq. (I2.9j) may be replaced 
by higher-dimensional operators such as q^d^H (0 / 'Mheavy) p ', where p is an integer and 
-^heavy is the mass scale where this dimension- (4 + p) operator is generated. Since the 
weak-singlet scalar <fi has a nonzero charge under U(l) z , the relations between the various 
quark charges may be changed. The higher- dimensional operators may be induced in a 
renormalizable quantum field theory by the exchange of heavy fermions that have Yukawa 
couplings to both the Higgs doublets and <p. Another loophole is that both Higgs doublets 
may be charged under U(l) z if there is a fine-tuning as discussed in section 2.1, so that the 
restrictions on quark charges would be again modified. Based on these considerations, 
we will study in some detail the implications of Eq. (|2.10j) . but we will also consider 
departures from it. 

We point out that generation dependent quark charges lead to flavor- changing cou- 
plings of the Z' in the mass eigenstate basis, where the fermion mass matrices are diagonal. 
Various experimental constraints from flavor- changing neutral current processes impose 
severe constraints on such flavor- changing Z' couplings, unless Mz> is so large that the 
effects of such a Z' would be beyond the reach of even the LHC To avoid these com- 
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plications we will avoid generation dependent quark charges in this paper. In practice 
this does not restrict significantly the generality of our results because the Tevatron is 
typically not very sensitive to Z' decaying into quarks, and the production cross section 
depends only on an average quark charge. Therefore, altogether there are three quark 
charges relevant in what follows: z u , Zd, z q . 

The masses of the electrically-charged leptons can be induced by the last term shown 
in Eq. ()2.9|) even if the lepton ^-charges are generation dependent. Moreover, no flavor- 
changing neutral currents are induced in the lepton sector by Z' exchange. The lepton 
mass terms impose, though, a relation between the left- and right-handed lepton z-charges: 
zn = z ej , j = 1,2,3, or zl — z ej = Zn x in the two-Higgs doublet model with large tan/3. 
As in the case of the quarks, we allow for deviations from these equalities motivated by 
lepton mass generation via higher-dimensional operators. Thus, all six lepton charges, 
zi., z ej , j = 1, 2, 3 could be relevant for Tevatron studies. 

Additional constraints arise due to the requirement of generating neutrino masses and 
mixings. The terms in the Lagrangian responsible for these are given by 

^Fi/^ T W2 /f + A^Ji^ 2 4'ift + m ^,^4' + h . c . , (2.11) 

where we have included right-handed neutrinos, z/^, which are singlets under the SM 
gauge group. If these are not present, then the last two terms in the above equation 
vanish. If there are n right-handed neutrino flavors, then j f , k' = 1, ...,n. For n > 2 all 
dimensionless coefficients Cjk of the above lepton-number-violating terms may vanish. The 
other parameters appearing in Eq. ()2.11j) are as follows: M v is the mass scale where the 
lepton-number-violating terms are generated, vn v -, k , are right-handed neutrino Majorana 
masses, AJ fc , are some Yukawa couplings. 

The requirement that the three active neutrinos mix, so that the observed neutrino 
oscillations can be accommodated, implies that the lepton charges are generation inde- 
pendent. However, as in the case of quarks, the terms in Eq. (j2.11|) may be replaced 
by higher- dimensional operators involving powers of 0/Mh eaV y Furthermore, the tiny 
neutrino masses make the existence of such higher- dimensional operators an attractive 
possibility [Sj. If several scalars carry different U{l) z charges, one could avoid almost 
entirely the constraints from neutrino mixing on lepton charges. 

The six lepton charges determine the leading decay width of the Z' into the corre- 
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sponding leptons: 



T(Z' 
T(Z' 

where ej = {e, /x, r} for j = 
decay widths of the Z'\ 

r(Z'->jets) 
F(Z' -> 66) 

r(z' -> «) 



' z f m z , 



(4 + 4+4) 



24tt 



-My, 



[2.12) 



1,2,3. Similarly, the quarks charges determine the following 



(24 + 4 + 4)f^, 



:4 + 4)lz^ fi + - 



[4 + 4)f Mz> 



7T 



Am z t 
Ml, 



2\l/2 



1 + ^ + (a s m 2 t /M 2 z ,)] 9 {M z , - 2m t ) 



(2.13) 



where "jets" refers to hadrons not containing bottom or top quarks and we have included 
the leading QCD corrections, but we have ignored electroweak corrections and all fermion 
masses with the exception of the top-quark mass, m t . Additional decay modes, into pairs 
of Higgs bosons (if z Hl ^ or zh 2 0), CP-even components of the scalar, right-handed 
neutrinos, or other new particles might be kinematically accessible and large. Therefore, 
the total decay width, IV, is larger than or equal to the sum of the seven decay widths 
shown in Eqs. (|2~P2J) and (|2~T3j) . 

Assuming that the decays into particles other than the SM fermions are either invisible 
or have negligible branching ratios, the Z' properties depend primarily on eleven param- 
eters: mass (Mz>), total width (IV), and nine fermion couplings (z ej , zn, z q , z u , zj) xg z . 



2.3 Realistic models 

So far we have imposed SU(3)c X SU(2) W x U(l)y x U(l) x gauge invariance on the La- 
grangian. Additional restrictions need to be imposed in order to preserve gauge invariance 
in the full quantum field theory: the fermion content of the theory has to be such that 
all gauge anomalies cancel. In our case, we need to make sure that there are no gauge 
anomalies due to triangle diagrams with gauge bosons as external lines. 
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Triangle diagrams involving two gluons or two SU{2) W gauge bosons, and one U{l) z 
gauge bosons give rise to the [SU(3)c] 2 U(l) z and [SU(2) w ) 2 U(l) z anomalies: 

A 3 3z = 3 (2z q - z u - z d ) , 

3 

A 2 2z = 9z q + J2 z h ■ ( 2 - 14 ) 

3=1 

Triangle diagrams involving U(l)y x U(l) z gauge bosons give rise to the [U(1)y] 2 U(1) z , 
U(1) y [U(1) z } 2 and [U(l) z ] 3 anomalies: 

3 

A Uz = 2z q - lQz u - Az d + 2 ( z h ~ 2 ^) » 

3=1 
3 

A Uz = Q(z 2 q -2z 2 u + z 2 d ) - 2^ (zl .-4.) , 

3 n 

A zzz = 9(2zl-z 3 u -z 3 d )+J2(H-4 3 )-J2i> ( 2 - 15 ) 

j=l i=l 

where we have included n right-handed neutrinos of charges z Vi under U(l) z . Finally, 
triangle diagrams involving two gravitons and one U(l) z gauge boson contribute to the 
mixed gravitational-C/(l) z anomaly, which makes general coordinate invariance incompat- 
ible with U(l) z gauge invariance: 

3 n 

A GGz = 9{2z q -z u -z d )+J2{ 2z ij- z e,)-J2 z ^ ( 2 ' 16 ) 

3=1 i=l 

Gauge invariance at quantum level requires that all the anomalies listed in Eqs. ()2.14|) - 
(I2.16|) vanish, or are exactly canceled by anomalies associated with some new fermions 
charged under both the SM gauge group and U(l) z . The impact of the new fermions on the 
Z' properties described here can be ignored if they are heavier than the Z'. Altogether, 
there are six equations that restrict the nine ^-charges of the SM fermions. Finding 
solutions to this set of equations is a nontrivial task, especially if one imposes that the 
charges are rational numbers, as suggested by grand unified theories. 

The case where the top-quark mass is generated by a Yukawa coupling to a Higgs 
doublet, as in the SM or the Minimal Supersymmetric Standard Model (MSSM), leads to 
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Eq. (|2.1U|) . in which case the A 33z , A 22z and Au z anomalies vanish only if 

1 3 

z q = z u = z d = -- 22 z h > ( 2 - 17 ) 

3 

Ek-^) = °- ( 2 - 18 ) 

The remaining anomaly cancellation conditions, in the absence of exotic fermions, take 
the following form: 

£(4-0 = »> 

3=1 

3 n 

EH -4) =£4 

i=i i=i 

(2.i9) 

i=l 

It is hard to find general solutions to this set of equations. A well known nontrivial 
solution is n = 3 and z\. = z ej = z v . = Z[, j = 1,2,3, which corresponds to the U{1)b-l 
gauge group. The associated gauge boson, Zb-l-, is an interesting case of a "non-exotic" 
Z' [H] relevant for Tevatron searches. We have found a generalization of this solution 
which preserves the Z\ . — z e . — z v . equalities within each generation, but have different 
lepton charges for different generations. In particular, the case of z^ = is worth special 
attention because the Z' does not couple to electrons, so there are no tight limits from 
LEP. In this case there are only two independent charges: z q which sets all quark charges, 
and Zi 2 which sets the charges of the muon and second-generation neutrinos. Normalizing 
the gauge coupling such that z q = 1/3, the r and third-generation neutrinos have charge 

z h = -3 - z l2 . (2.20) 

If new fermions are included, the anomaly cancellation conditions have more solutions. 
We have found that including within each generation two fermions, ip l and ip e , which under 
the SM gauge group are vector-like and have the same charges as II and e#, respectively, 
allows z charges proportional to B — xL with x arbitrary. The U(1)b- x l charges are 
shown in Table 1. This is the most general generation-independent charge assignment 
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Table 1: Fermion gauge charges. 



for the SM fermions that allows quark and lepton masses from Yukawa couplings, and is 
relevant for Z' searches at the Tevatron. 

All other generation-independent U(l) z charge assignments require the restrictions 
on fermion charges from fermion mass generation to be lifted, for example by replacing 
the Yukawa couplings with higher- dimensional operators. The six anomalies given in 
Eqs. (I2.14|) - (J2.16|) vanish only for the nonexotic family of U(l) z charges that depends on 
two parameters pj- Assuming that z q ^ 0, and normalizing the g z gauge coupling such 
that z q = 1/3, determines all other charges as a linear function of a single free parameter, 
x, as shown in Table 1. We label this charge assignment by U(l) q+xu . Particular cases 
of Z' "models" include U(1)b-l for x = 1, the U(l) x from SO (10) grand unification for 
x = —1, and the [U(1)r x U(1)b-l]/U(1)y group from left-right symmetric models for 
x = 4 — 3g 2 R / g Y where g R is the U(1)r gauge coupling. 

Many popular Z' models are accessible at the Tevatron provided both the restrictions 
from fermion mass generation are lifted and new fermions charged under the SM gauge 
group are present. We have found a couple of generation- independent charge assignments 
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that depend on a free parameter x, which include the ^-inspired Z' models that have been 
frequently analyzed by experimental collaborations. Both require within each generation 
two fermions, if and ip d , which under the SM gauge group are vector-like and have the 
same charges as II and cIr, respectively. One assignment, labeled by U(l)d- xu , has z q = 0, 
z d = 1/3 and z u = —x/3. For x = the i? 6 -inspired U(l)j is recovered, while x = 1 
gives the "right-handed" U(1)r group. The other assignment is such that all fermions 
belonging to the 10 representation of the SU(5) grand unified group have the same U(l) z 
charge, assumed to be nonzero and normalized to 1/3, while the fermions belonging to 
the 5 representation have charge x/3. We label this assignment by U(l) w+X ^. Anomaly 
cancellation requires two right-handed neutrinos per generation. Particular i? 6 -inspired 
cases include U{1)^ for x — 1, U{l) x for x = —3, and 11(1)^ for x = —1/2. Note that when 
the LEP and Tevatron experimental collaborations refer to these particular models, the 
gauge coupling is usually assumed to be determined by a unification relation: g 2 z = gy£, 
with f = 5/8 for U{l) f , £ = 3/8 for U{l) x , £ = 1 for U(l) n , £ = 5/3 for U{l)i. Thus, our 
families of models completely describe the physics of the GUT- inspired £7(1) 's, but also 
allow one to relax their assumptions and explore more general Z' physics. 

There is an interesting class of Z' models in which the Higgs is a pseudo-Goldstone 
boson of a spontaneously broken global symmetry. These "Little Higgs" models always 
include at least one Z' to cancel the leading quadratic divergence in the Higgs mass from 
loops of the ordinary W and Z. A proto-typical model of this type is the "Littlest Higgs" 
|12| . This already reveals a key feature of the little Higgs Z'\ it always couples to the 
Higgs, and thus generically has strong constraints from Z-Z' mixing, requiring the Z' 
mass to be larger than several hundred GeV [13]. Thus, it is usually not very interesting 
for Tevatron searches. 

2.4 LEP II limits on Z' models 

The constraints from e + e~ colliders on the Z' properties fall into two categories: precision 
measurements at the Z pole, through the Z-Z' mixing discussed previously, and measure- 
ments of e + e~ — > ff above the Z-pole at LEP-II, where / are various SM fermions. In 
practice, the agreement with the SM requires that either the Z' gauge coupling is smaller 
than or of order 10 -2 jS], or else Mz* is larger than the largest collider energy of LEP-II, 
of about 209 GeV. In the latter case, of interest at the Tevatron, one can perform an 
expansion in s/M|,, where s is the square of the center-of-mass energy. This leads to 
effective contact interactions which have been bounded by LEP-II for all possible chiral 
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structures and for various combinations of fermions. These interactions are parameterized 
by the LEP electroweak working group [Hj as, 

( 1 + M(A^ W (/7 " Pb/) (2 ' 21) 

where Pa is a projector for right- (A = R) or left-handed (A = L) chiral fields, and 
8 e f = 1, for / = e, f 7^ e, respectively. These contact interactions provide a model- 
independent framework in which LEP-II data can confront high mass effects beyond the 
SM, up to corrections of order s/M|,. 

In the absence of flavor violation in the Z' couplings, the Z' contributions to e + e~ — > 
// for / ^ e proceed through an s-channel Z' exchange, with tree level matrix element, 

jj^- s [eiMPL + z e P R )e] [h^z fL P L + z fR P R )f] (2.22) 

where tiny terms of order m/m e /M§, have been dropped. In the case of / = e, there is 



also a t-channel exchange, which in fact motivates the factor of (1 + 5 e f) in Eq. (j2.21j) . 
which allows one to treat all / equivalently at the level of matrix elements. 

One should compare the LEP A LL , A RR , A LR , and A RL limits to the operators of each 
structure in the Z' theory in order to find a limit on a given Z' model. This procedure 
finds the best single bound from each operator on a given Z' model, but it ignores the 
potentially stronger bound that comes from the combined effect of more than one operator. 
In the absence of a dedicated reanalysis of the data, this is the best one may do. However, 
we reiterate that it does not always represent the best potential bound from the data, 
due to correlated effects on observables which cannot be taken into account correctly in 
this way. Typically the strongest bound comes from a single choice of chiral interaction 
combination and /. 

Matching the Z' matrix elements, Eq. ()2.22)1 . to the LEP-II formalism, Eq. (j2.21j) . one 
derives bounds such as, 

M\,-s> ^\z eA z fB \{A^ B )\ (2.23) 

and one must choose A + for z eA zj B > or A~ for z eA zj B < 0. Typically, the LEP-II 
bounds on the A's are on the order of 10 TeV, schematically translating into bounds on 
the Z' mass on the order of Mz> ^ \z\g z x (a few TeV). More precise bounds in the 
case of certain models are discussed below. The Tevatron can effectively improve our 
knowledge of Z' models only when the couplings zg z are appropriately small such that 
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Figure 1: Lower bounds on Mz>/g z from the LEP-II search for LL, RR, LR and RL 
contact interactions, applied to the models of Table ^ as a function of the continuous 
parameter x. For U(1)b-xL, we have included the bound on vector-like e + e~ — > £ + £~, as 
is appropriate for that model. 

the LEP-II limits are evaded for Z' masses in the several hundred GeV range, but large 
enough that the Z' rate is observable compared to backgrounds. 

For example, U(1)b- x l has vector-like interactions with quarks and leptons, and thus 
it is better to compare LEP limits for vector-like interactions (Ayy) than for the individual 
chiral set of All, Arr, Alr, and Arl- From Ref. [H], we see that the strongest of the 
Ay V bounds is from the process e + e~ — > £ + £~: A + > 21.7 TeV. This is in fact one of the 
most stringent bounds to be found from LEP-II. Translating this bound in the specific 
couplings of U(1)b-xL results in the limit Mz> > \x\g z x (6 TeV). For the remaining 
models of Table ^ the analysis is more complicated, because the best bounded channel 
typically depends on the value of x. Thus, for these models, we have scanned (for fixed 
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—10 < x < 10) through all channels of A LL , A RR , A LR , A RL and chosen the best bound. 
The results are shown in Figure [T] It is interesting to note that for \x\ > 1, U{1)b- x l is 
very strongly constrained by LEP II data, whereas for x — > the coupling to electrons 
becomes small and the bounds disappear. 

We have compared our results for general x with the ^-inspired models studied in 
detail at LEP-II [Hj. As explained above, these correspond to specific points in x for a 
given model family. We find that our results agree with the LEP bounds for the dedicated 
analysis at or better than the 25% level (depending on the model), thus indicating that 
our procedure does a good job in comparison with the dedicated analysis for the points in 
which the two may be compared. For most of the parameter space there is no dedicated 
LEP analysis, and we present for the first time the LEP bounds on the general class of 
Z' models. 

Generically, the fact that LEP II was an e + e~ collider implies that these strong bounds 
can be evaded by a Z' which couples only very weakly to electrons 2 . Also, different chiral 
structures than the vector-like interactions of U(1)b- x l can have weaker bounds. From 
Ref. JH] one observes that a Z' which couples only to left-handed or right-handed electrons 
is bounded only by A + > 7.1,7.0 TeV (for All and A RR , respectively). This implies a 
bound of Mz> > g z z x (1-9 TeV), approximately three times weaker than the bound on 
U(l) B - xL . 

3 Z' Searches at the Tevatron 

At the Tevatron, searches for additional neutral gauge bosons can be performed in a 
variety of processes. If such bosons couple to the SM quarks, they may be directly 
observed through their production and subsequent decay into high energy lepton pairs or 
jets. The case of the decay into leptons is particularly attractive due to low backgrounds 
and good momentum resolution. Bounds on several models containing extra neutral gauge 
bosons, have been set by both the CDF [TSJ EH 1201 and DO |2H 122 123 experiments by 
measuring high energy lepton pair production cross sections. Searches have been made in 
the e + e~ channel, which has the best acceptance, and thus best systematics, as well as in 
the channel. More recently, the challenging t + t~ channel has also been analyzed 

|24j . The /i + yU~ and t + t~ final states, along with the Z' decay into jets which suffers from 

2 For example, theories such as Top-color assisted Technicolor Top-flavor |16| . or Supersymmetric 
Top-flavor have small couplings to the first generation. 
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huge QCD backgrounds, can probe Z 1 bosons with suppressed couplings to the electrons, 
which are not constrained by the LEP searches. 

In what follows, we will restrict ourselves to the study of the leptonic decay modes, 
proposing a simple, model- independent, parameterization for the Z 1 production cross 
section and analyzing its theoretical and experimental feasibilities and limitations. 

3.1 Z' Hadro-production 

The additional terms, beyond those coming from SM particles, in the differential cross 
section for production of a pair of charged leptons due to the presence of an extra neutral 
gauge boson can be written as 25 j 

-^a (pp -> Z'X - l+l-X) = - s o{Z> - l+r) W Z , (s, Q 2 ) + ^ , (3.1) 

where Q is the invariant mass of the lepton pair, and y^s is the energy of the pp collision 
in the center-of-momentum frame. The first term accounts for the contributions from 
the Z' itself and has been explicitly factorized into a hadronic structure function, Wz>, 
containing all the QCD dependence and the couplings of quarks to the Z', and 

gl (zl+Z 2 e \ Q2 

a{Z ' _ H) = -L \^^- j ( g 2 _ M 2 ;)2 + M | ;r | ; • ( 3 - 2 ) 

Up to NLO in QCD, only the partonic processes qq — > Z'X (non-singlet) and qg — > Z'X 
contribute to the hadronic structure function. If the Z' couplings to quarks are generation 
independent, both processes give contributions which are proportional to (z 2 + z 2 ) or 
( z q + z d) f° r U P anc ^ down type quarks, respectively. Therefore, the hadronic structure 
function can be written as 



W z > {s,Ml,)=gl 



{z\ + z 2 u ) w u (s, M 2 Z ,) + (z 2 q + z\) w d (s, M 2 Z ,) 



(3.3) 



The functions w u and Wd do not depend on any coupling and are exactly the same for 
any model containing neutral gauge bosons coupled in a generation independent way to 
quarks. In the MS scheme, they are given by 

I dxi I dx2 I dz{f q/P {x 1 ,M 2 z ,)fg /P {x 2 ,M 2 z ,)A qq {z,M 2 Zf ) 



q=u,c(d,s,b) 



+ f g/P (x u M 2 ,) [f q/P (x 2 , M 2 ,) + f q/P (x 2 , M 2 ,)] A gq (z, M 2 ,) 



Ml, 



+ ( Xl ^x 2 ,P^ P)\5 l^-zx 1 x 2 , (3.4) 
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where fi/p(x, M 2 ) and fi/p(x, M 2 ) are the PDFs for the proton and antiproton, respec- 
tively, and 



A qq (z,M*,) = 5(l-z) + ^^-C F 



71 



l+z 2 

2 (1 + z) ln(l -z)- In z 



2tt 



(1 - 2z + 2z 2 ) In. 



+ - + 3z - -z 2 



-*[o.l] 
(3.5) 

(3.6) 



The color factors are Cp = 4/3 and Tp = 1/2, and we have set the renormalization and 
factorization scales to Mz>- 

The second term in Eq. (|3.1|) . da- m i/dQ 2 , corresponds to the interference of the Z' 
with the Z and the photon. If the Z' resonance is narrow enough, the interference of the 
Z' with the Z and photons can be neglected (see the Appendix). 

In the narrow width approximation, the expression for the total cross section is simply 
obtained from the differential cross section, explicitly 

TT 



a (pp -> z'x -> i+rx) 



48 s 



W z > (s, M|,) Bt(Z' -> i + r 



(3.7) 



where By(Z' — > / + / _ ) is the branching ratio for the decay of Z' into the corresponding 
pair of leptons. Using the expression of the hadronic structure function, Eq. ()3.3j) . one 
obtains 



7T 

48s 



[c u w u (s,M|,) +c d w,i (s,M|,)] 



The coefficients c u and q, given by 
c u4 = g 2 z (z 2 + z 2 u4 ) Br(Z' - Z+i" 



(3.8) 



(3.9) 



contain all the dependence on the couplings of quarks and leptons to the Z\ while w u and 
Wd only depend on the mass of the gauge boson and can be calculated in a completely 
model-independent way. 

The parameterization given in Eq. ()3.8|) permits a direct extraction of a bound in 
the c u — Cd plane from the experimental limit for the cross section, which can be later 
compared to the predictions of particular models. This fact is particularly useful for 
models admitting free parameters like the ones discussed in the preceding sections. In 
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Figure 2: Excluded regions in the q — c u plane from the current 95% CL. limit for 
a ■ Br(Z' — > Z + Z _ ) given in [2H], for different values of the Z' mass. The thick straight 
line corresponds to values of c u and q in the B — xL model, in which c u = Cd- The area 
between the two thin straight lines is the region where the q + xu model lies. 

particular, these quantities are simply computed for a given Z' model, without need to 
compute the hadro-production cross section, and thus are a common ground between 
theory and experiment. 

The DO and CDF Collaborations have set preliminary 95% C.L. limits for a ■ Br(Z' — > 
in Run II with 200 fb _1 |27| |2H]- in Figure 121 we show the excluded regions in the 
Q — c u plane for different values of the mass of the Z' boson as obtained from the limit 
on a ■ Br(Z' -> l+l~) given by the CDF Collaboration [2E] in Run II with 200 fb _1 . Very 
similar results are obtained using the results by the DO Collaboration [2Zj. The w u and 
w d coefficients in Eq. flUBl) for this plot were calculated at NLO with MRST02 PDFs [2H]- 
From the current generation of CTEQ 30 PDFs, one obtains very similar results. 

It is instructive to compare these limits with the predictions of the four families of 
models presented in Section 2. The values of c u and q as functions of the gauge coupling 
g z and the x parameter are given in Table 2 In Figure El displays the values of (q, c u ) 
corresponding to the B — xL and q + xu models. In the B — xL case, these points are 
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1+x 2 
2 


J_ 

X 2 



Table 2: Predictions for a and c u in four families of models defined in Table 1. The 
branching fractions are computed at tree level for M' z > 2m t and assuming decays 
only into SM particles. 

constrained to satisfy c u = q, corresponding to the thick straight line. For the q + xu 
model, the allowed region is 

(3 - 2V2) C(i < c u < (3 + 2^2) c d , (3.10) 

which corresponds to the area between the two thin straight lines in Figure|21 The 10 + 2 5 
model, in turn, is constrained to the region c u < 2 q, whereas there are no constraints for 
the possible values of q and c u in the d — xu model. 

3.2 Higher-order corrections 

At NNLO in QCD, the general expression in Eq. (|3.3|) is no longer valid. This is due to 
the contributions from the partonic processes qq — > Z'X (singlet) and qq — > Z'X, which 
depend upon a variety of coupling combinations in addition to [z 2 + z 2 ) and (z 2 + z 2 d \ 
Thus, one should worry whether c u and q are a sufficient description of the model. The 
actual size of these corrections can be estimated by looking at a particular model. In 
Figure El we plot the sizes of the 0(a s ) and 0(a 2 ) terms to the structure function at 
NNLO relative to the Born contribution for the case of SM-like couplings of the Z' boson. 
The NNLO corrections were calculated with the program ZPROD [2*51 l2*o1 13T] and the 
MRST02 NNLO set of PDFs [H]. We have split the 0(a 2 ) corrections in two parts, 
one proportional to q and c u , and the other depending upon other combinations of the 
couplings. Contributions proportional to q and c u , coming from 0(a 2 ) corrections to 
processes already present at lower orders, are clearly the dominant ones, overcoming the 
remaining pieces by more than an order of magnitude in the whole Q range. Typically 
the terms with mixed couplings contribute less than one per mil to the structure function, 
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Figure 3: The ratio of the 0(a s ) and 0(a 2 ) corrections to the hadronic structure function 
over the Born contribution, assuming SM couplings for the Z' . The full line corresponds 
to the 0(a s ) corrections. The dashed and dotted lines, which, within the resolution of 
the figure, appear as a single line, correspond to the total 0(a 2 ) contributions and to the 
0(a 2 ) corrections retaining only those terms that give contributions proportional to c u 
and Cd to the cross section, respectively. 

while the total 0(a 2 ) corrections amount between 2% and 20% [for comparison, the 0(a s ) 
terms contribute between 20% and 50% of the structure function]. 

Although the actual values of the different higher order corrections will depend upon 
the model considered, it is reasonable to expect that terms not contributing to the Q-c u 
piece of the cross section will be negligible at NNLO. In that case, the parameterization 
in Eq. (J3.3)) holds as a very good approximation, and experimental bounds can be set in a 
model independent way, taking q, c u and Mz> as the only relevant parameters at NNLO 
accuracy. 

It is customary to take into account higher order corrections to the structure func- 
tion by means of a fc-factor, allowing the use of LO Montecarlo simulations, corrected 
afterwards by the mentioned factors. As the structure function at LO has a different 
dependence upon Q than at NLO or NNLO, the corresponding fc-factors, defined as the 
ratio of the higher order results over the LO one, k N i LO = W^ LO /W^P , vary noticeably 
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Figure 4: NLO and NNLO /c-factors for SM-like couplings as a function of the invariant 
mass of the lepton pair. Also shown is an approximation, Eq. (|3.11j) . for the NLO /c-factor 
which only takes into account the soft corrections to the structure function in the DIS 
scheme. 

with the invariant mass of the lepton pair. This variation has to be properly taken into 
account when correcting LO Montecarlo results. 

In Figure |U we plot the NLO and NNLO /c-factors obtained again with the program 
ZPROD, for SM-like couplings. There, we also show a usual approximation for the NLO 
/c-factor (used, for instance, in Ref. [20]). which takes into account only the soft pieces of 
the NLO structure function in the DIS scheme, namely 



The variation of the /c-factor, in the plotted range of Q, for the NLO case amounts to more 
than 10%, the genuinely NNLO corrections are also sizable whereas the soft approximation 
does not provide a good description of the higher order effects. The difference between 
the soft approximation and the full /c-factors, yields a discrepancy of about 5% to 10% 
for the cross section in the high mass region, which is particularly relevant in the search 
of extra gauge bosons. Similar results are obtained using a flat /c-factor with k = 1.3, as 
in reference [28 , instead of the soft approximation. 




(3.11) 
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So far we only discussed QCD corrections to the Z' cross section, which are of the 
order of 30% as can be seen from the fc-factors in Figure EJ There are also corrections from 
the electroweak sector, which we will address briefly. The complete 0{a) corrections to 
the SM contributions to neutral current Drell-Yan process were calculated in (321 . There, 
it was found that these corrections are large, particularly in the high invariant mass 
region, being of the order of 12% at Tevatron energies for ran — 700 GeV in the electron 
channel. Besides affecting the background for the search of additional neutral gauge 
bosons, electroweak radiative corrections also modify the signal cross section. However, 
as we will see, these effects are substantially smaller for the Z' terms. 

As shown in [32] , the main contributions to the electroweak corrections come from the 
box diagrams and cannot be factorized into effective couplings and masses. In particular, 
box diagrams with two charged bosons give rise to large double logarithms which are 
the origin of the large corrections in the high mass region. On the other hand, box 
diagrams that include neutral bosons (7, Z° and Z') always appear in combination with 
their crossed versions and that leads to a cancellation of the double logs [33]. Then, the 
non-factorizable contributions that affect exclusively the signal cross section only include 
subleading simple logs and thus the corrections are expected to be smaller than for the 
SM background. 

The remaining contributions come from QED and factorizable purely weak corrections. 
The later can always be absorbed into effective couplings and masses, thus, they do not 
affect the signal cross section where these quantities are treated as free parameters. The 
main electromagnetic contributions come from large logarithms due to collinear photon 
emission in the initial and final states, and affect both the signal and background cross 
sections. The large contributions coming from initial state radiation can be factorized 
into the PDFs, modifying the DGLAP evolution equations for the partonic densities. 
After factorization, the remaining terms are typically at the per mil level, reaching 1% in 
the high momentum fraction region [H2] whereas the QED modifications to the evolution 
equations are small and neglected in comparison to the uncertainties in the PDFs [8lj . 
However, collinear emission in the final state gives corrections of the order of 5% for the 
electron channel in the high mass region, and a careful analysis should probably take 
them into account. 
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3.3 Model dependence in experimental bounds 

As we have shown in the previous section, the parametrization given in Eq. (|3.8|) allows to 
extract model independent constraints on the coefficients q and c u from the experimental 
results for the lepton pair production cross section. A key assumption for this analysis is 
that the bounds for the cross section can be extracted from data in a model independent 
manner. In particular, the experimental analyses involve corrections for the finite accep- 
tance of the detectors to extract the total cross section. As the acceptance is obtained 
from detailed Monte Carlo simulations, which need to assume particular values of the 
couplings, it is far from trivial that this procedure does not introduce model dependence 
into the experimental bounds. 

In Ref. [TBI, ^he changes in the acceptance with variations in the couplings of up and 
down quarks were studied. There, it was found that the acceptance changes very little 
when considering the limiting cases of either decoupling up or down quarks. However, 
this study was limited to SM-like couplings for electrons, a feature that, a priori, might be 
too restrictive. To study the actual model dependence of the experimental acceptance, in 
this section, we will consider the angular distribution of the lepton pair and apply simple 
cuts on this distribution. For simplicity we will restrict to the LO approximation. 

In the left panel of FigureElwe plot the LO cross section differential with respect to the 
azimuthal angle, in the center of mass frame of the lepton pair, with different assumptions 
for the couplings, setting Mz> = 600 GeV. We considered SM-like couplings, SM-like 
couplings with up or down couplings neglected and the E 6 inspired model U(l)j mentioned 
in Section 2. The right panel shows the ratios between the cross section in the last three 
cases to the cross section with SM couplings. Except for the overall normalization, the 
cases where the Z' does not couple to either up or down type quarks differ very little from 
the SM case. This feature can be traced back to the peculiar fact that, in the SM, the 
left and right handed couplings of charged leptons, satisfy zf — z\ = sin 2 9w — 1/4 — 0.02. 
Then, the terms odd under 9 — > —9 are suppressed relative to the even ones in the LO 
differential cross section, which turns out to be nearly symmetric. This characteristic 
feature is in sharp contrast with the behavior in the U(l)j case, where the asymmetry 
is almost maximal. This is related to the vanishing of the couplings of the right-handed 
electrons and left-handed down-type quarks to the Z' in this last model. 

To get a handle on how the noticeable differences in the angular distribution affect 
the experimental acceptance, we crudely estimated it by integrating the differential cross 
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Figure 5: (a) Angular distribution, in the CM frame of the lepton pair, at LO in QCD 
for different models. The lines labeled no u (d) coupling correspond to SM charges but 
neglecting the coupling of up (down) type quarks to the Z'. (b) Ratio between the 
differential cross section (shown in panel a) in different models and the case of SM-like 
charges. 

section imposing a symmetric cut on the lepton angle, defined in the laboratory frame, 
and normalizing to the total cross section. In Figure |H1 we plot the results obtained 
for the models considered in the previous paragraph, taking the U(l)j case as reference 
normalization. Although the angular distribution in this model differs substantially from 
that in the the other cases, the corresponding acceptances coincide at the few percent level. 
Note that the SM case without down-type quark couplings is practically indistinguishable 
from the U(l)j one. 

This, apparently peculiar, result for the acceptance is due to the fact that, on average, 
the vector boson is produced with very small longitudinal momentum, so the boost to 
go from the laboratory frame to the center of mass one is small. Thus, a symmetric cut 
in the lab frame corresponds to an almost symmetric cut in the center of mass frame, 
which in turn means that the ratio we are using to estimate the acceptance has only 
small contributions that depend on the couplings of the quarks and leptons to the Z' . 
In particular, models that do not couple at all to up (down) type quarks give practically 
identical results for the acceptance. 

In conclusion, we find that, even though the angular distribution of the final state 
leptons is highly model dependent, the experimental acceptance, which we naively esti- 
mate through angular cuts, is only mildly affected by this dependence, in accordance with 
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Figure 6: The Acceptance (A), computed from the LO cross section integrated over the 
azimuthal angle of the leptons in the laboratory frame in the region 50° < 0^ < 130°. 
We considered the cases of SM couplings, SM couplings without w-quark couplings and 
without ci-quark couplings respectively. To show the variation compared to a most extreme 
situation, we normalize to the corresponding LO acceptance of the U(l)i model. 

previous studies by the experimental groups. On the other hand, the model dependence 
of the angular distribution potentially can be a tool to discriminate between models in 
case of a discovery. 

3.4 Probes of Z' models at the Tevatron 

With the forthcoming data, the Tevatron experiments will be able to explore a new region 
in masses and couplings for many models containing extra neutral bosons, with the chance 
of making a Z' discovery. 

In case that no signal excess is observed, the extracted limits to the Z' cross section 
will set new bounds on masses and couplings, or alternatively on the c u and Cd parameters 
discussed in previous sections. As an illustration, in Figured we show estimates for the 
experimental limits, for different values of the integrated luminosity. For this estimates, 
we considered the current CDF bound [2H] and assumed that the limit scales as the inverse 
of the square root of luminosity, as will be the case provided the limit is dominated by 
statistics and not systematics. In the two upper panels, we plot the bounds on cross 
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Figure 7: Projected bounds for the B — xL (left upper panel) and 10 + x5 (right upper 
panel) models and projected excluded regions in the c<j — c u plane (lower panel). In the two 
upper panels, the vertical marks show the current LEP bounds for the Z' mass obtained as 
described in Section 2.3. In the B — xL case, for x = 3 this last bound is Mz> > 1800 GeV; 
for the 10 + x5, the bound for x = is Mz> > 119 GeV, beyond the scope of the figure. 
The projected bounds at luminosities C = 2fb _1 , 10fb _1 are obtained from Fig. 2 by 
scaling with a factor 1 / \JH. The lower panel also shows the regions in the Cd — c u plane 
corresponding to the B — xL and q + xu models, as in Figure El showing the projected 
increase in reach with 2 fb~ x . The dots labeled a, b and c correspond to the B — L model 
with g z = 0.1, g z = 0.3 and g z = 0.5, respectively. 

section times branching ratio, as a function of Mz>, together with the predictions for 
different values of x in the B — xL (left) and in the 10 + x5 (right) models. We also show 
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the mass bounds, for the different cases, set by the LEP contact interaction constraints 
discussed before. The lower panel shows the excluded regions in the q — c u plane for 
different values of the Z' mass, assuming an integrated luminosity of 2fb~ 1 . 

The results in Figure[7|show that, both in the cases of the B — xL and 10 + x5 models, 
there is a sizable unexplored region in parameter space that the Tevatron will certainly 
be able to probe. For the B — xL models, LEP bounds are stronger for larger values of 
\x\, while the Tevatron can do much better for \x\ < 1. For the 10 + a;5 models, the LEP 
bounds are slightly weaker than in the previous case and the unbounded region in x space 
is larger. 

If a signal excess is seen in the invariant mass distribution, it would also be possible to 
shed some light on the nature of the couplings of the new boson to the fermions by studying 
the angular distribution of the final state leptons. As shown in the previous section, there 
are substantial differences, between models, in the predicted angular distribution. 



CO 

a- 



-0.5 



M z ,=700GeV 



• CDF data 
— SM 

°- 5 U(l) B . xL 
U(l) 




10 



M ee (GeV) 



Figure 8: The forward-backward asymmetry at LO for the B — L and U(l)d- X u models 
with M z > = 700 GeV. In the case of B — L we chose g z = 0.1. For the ^7(1)^-^ models, 
we fixed x = 1 and g z = 0.5 We also show the SM prediction and recent CDF data for 
this observable |35j . 

In Figure |H1 we show the forward backward asymmetry at the Tevatron, predicted in 
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the B — xL and d — xu models for the case of a Z' boson with Mz> = 700 GeV. The 
shift from the SM prediction of the forward backward asymmetry in these two models 
has opposite sign, allowing, in principle, to distinguish between them, provided enough 
statistics are collected in the high mass region. Observables at the Tevatron other than 
the forward-backward asymmetry are discussed in Ref . jSE] , while capabilities of the LHC 
and a high energy e + e~ linear collider are addressed in Ref jSj. 

4 Conclusions 

At the Tevatron, the hypothetical Z' bosons may be produced via their couplings to light 
quarks, and are more likely to be detected if they decay into charged leptons. The typical 
signature of a Z' boson would be a bump in the total cross section for dilepton production 
as a function of the dilepton invariant mass. The observability of a dilepton signal in the 
inclusive process pp — > l + l~X, where I stands for e, fj, or r, is controlled primarily by two 
quantities: the Z 1 mass, and the Z' decay branching fraction into times the hadronic 
structure function Wz> defined in Eq. ()3.3|) . The exclusion limits presented by the DO 
[21]- [22] and CDF Collaborations [111-120] are curves in the plane spanned by these two 
parameters. Such an exclusion plot is very useful, allowing one to derive the range of 
Z' parameters consistent with the experiment. However, the hadronic structure function 
entangles the model dependence contained in the quark- .Z 7 couplings with the information 
about the proton and antiproton structures contained in the PDFs. In order to simplify 
the derivation of the exclusion limits in the large Z' parameter space, we are advocating 
the presentation of the exclusion curve (see Figure 3) in the c u — q plane, where c u and 
Cd are the decay branching fraction into leptons times the average square coupling to up 
and down quarks, respectively. 

Assuming that the couplings of Z' to quarks and leptons are independent of the fermion 
generation, the Z' properties are described primarily by seven parameters: mass (Mz>), 
total width (T Z '), and five fermion couplings (z e , zi, z g , z u , xg z . The exclusion curve 
in the c u — Cd plane sets a bound on a single combination of these seven parameters. 
Nevertheless, in any specific model defined by certain fermion charges it is straightforward 
to compute c u and q, and to derive what is the limit on the gauge coupling g z as a function 
of M z >. 

For example, if the quark and lepton masses are generated by Yukawa couplings to 
one or two Higgs doublets, as in the SM or its supersymmetric versions, the only gauge 



27 



groups that may provide a Z' gauge boson accessible at the Tevatron are of the type 
U(1)b-xL- This means that all fermion charges are determined by a single parameter, x. 
Within this family of gauge groups, c u and q have a simple dependence on x and g z ; for 
a given x and M z >, the limit on g z can be immediately derived. 

If the quark and lepton masses are generated by a more general mechanism, Z' gauge 
bosons associated with gauge groups other than U(1)b- x l m ay be accessible at the Teva- 
tron. We have presented three other examples of one-parameter families of U(l) gauge 
groups, chosen to include (for particular values of the parameter x) many of the Z' models 
discussed in the literature. For these families of models, the Tevatron reach goes signif- 
icantly beyond the LEP II bounds for large regions of the three dimensional parameter 
space spanned by Mz>, g z and x. 

Relaxing the assumption that the couplings of Z' to leptons are generation inde- 
pendent, for each of the e + e~, fi + fi~ and t + t~ final states there is a different c u and 
Cd- Interestingly, U(l) gauge groups that lead to a Z' of this type exist even when the 
anomalies cancel without need for new fermions charged under the SM group, and the 
quark and lepton masses are generated by Yukawa couplings to a single Higgs doublet. 
Such Z' bosons may have very small couplings to electrons, evading altogether the LEP 
bounds, and could be discovered in the fi + fi~ or t + t~ channels at the Tevatron. 

Although generation-independent Z' couplings to quarks are tightly constrained by 
measurements of various flavor-changing neutral currents, a Z' with different couplings to 
the d and s quarks (or to the u and c quarks) in the mass range accessible at the Tevatron 
cannot be completely ruled out. In that case, the c u and q parameterization would have 
to be supplemented by c s (or c c , q,) quantities. Even in this case, the fact that current- 
and next-generation hadron colliders collide nucleons (and anti-nucleons) implies that c u 
and Cd typically remain the most important, because of their large valence distribution 
functions (particularly at large parton x) for nucleons. 

Observables other than the total cross section for dilepton production can also be 
measured at the Tevatron. We have discussed the additional information provided by 
the forward-backward asymmetry. In most cases however, a Z' discovery is more likely 
to occur first as a bump in the dilepton total cross section. If that happens, M z > can 
be determined by the invariant mass of the lepton pair, and a curve (actually a band of 
experimental error bars) in the c u — Cd plane can be derived. For each Z' model (fixed 
x within the one-parameter families), the curve would determine the gauge coupling. 
However, pinning down the model would be difficult, requiring additional observables at 
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the Tevatron and future colliders. 
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Appendix: Interference terms 

The Z' interference with the Z and the photon is taken into account by the second term 
in Eq. (JO}. This term can be factorized similarly to the contribution due solely to the 
Z' [first term in Eq. ()3.1j . by replacing a[Z' — > with 

v ,_9z9x ( z ljZ x + z ej z? \ (Q 2 - M 2 ,) (Q 2 - M\) + M z > M x Y z , Y x 

(Q 2 - M x f + M x T^j ' 

where X = 7, Z° and gx, {2*,%*), M x and T x are the corresponding coupling, lepton 
charges and mass and width of the boson. The quark charges in Wz> must also be changed 
accordingly in Eq. ()3.1|) . 

For a narrow Z' resonance, the interference of the Z' with the Z and photons can be 
neglected. As an illustration, in Figure^ we plot the NLO cross section for the production 
of an electron-positron pair as a function of the lepton system invariant mass. The curves 
shown correspond to the SM background for the process and to the Z 1 mediated ones, 
with and without the interference terms with the Z Q and the photon. For the plot we 
chose the B — xL model, discussed in the previous section, and fixed M z > = 600 GeV, and 
x = — 1. Note that for x = +1 the only difference is that the interference terms change 
sign. The two sets of curves correspond to different choices for the gauge coupling, namely 
g z = 0.05 and g z = 0.2. These values of the coupling correspond to T Z ' = 0.26 GeV and 
T Z t = 4.15 GeV, respectively, assuming that only decays to SM particles are allowed and 
neglecting all QCD and electroweak corrections. The bounds set by LEP for these two 
cases are M Z i > 300, 1200 GeV respectively. In the case of g z = 0.05, the signal cross 
section outside the Z' peak is completely negligible compared to the SM background, 
and, at the peak, the interference terms can be neglected. For g z = 0.2, the interference 
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Figure 9: NLO differential cross sections for the production of electron-positron pairs as 
a function of the invariant mass of the pair, in the U(1)b+l model. The solid curves 
correspond to the total cross sections for the signal, including interference terms, for the 
gauge coupling fixed to g z = 0.05 and g z = 0.2 respectively. Dashed lines are the same 
cross sections neglecting interference terms and the dotted lines correspond to the SM 
background. 

terms are more important and contribute to the tails at low and high mass. However, the 
experimental errors would not allow one to disentangle the signal from the background 
outside the peak, where again, the signal cross section is dominated by terms containing 
only Z' propagators. 
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